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Abstract 

We propose a new method of gravitational waves detection in the 
>S ! KT 1 -f- 1(T 2 - Hz band for a space laboratory based on the use of 

| the Kozorez effect in the magnetic interaction of superconducting 

O^ ■ solenoids. 
O 

cr 

£*'• 1 Introduction 

op. 

As Thorne noted QXj it can be expected that the amplitude h of gravitational 
^ ■ wave bursts from astrophysical sources reaches the values of the order of 

1CT 16 -i- 1CT 17 in 10—1 -T- 10~ 2 -Hz band . Some authors proposed to use 
in this frequency band the Doppler tracking of interplanetary spacecraft , 
Skyhook-detector and the exitations of seismic motions in the Earth's sur- 
face . A new type of the gravitational wave detector for the low frequency 
band is cosidered in this paper. The detector is based on the effect that 
the potential energy of a pair of magnetically interacting superconducting 
solenoids in weightless state , generally speaking, has the minimum at some 
nonzero distance. At this distance the solenoids are in a weak equilibrium 
condition. This system in vacuum is a weakly coupled nonlinear oscillator 
with very low energy dissipation and it can be used as a sensitive detector 
of tidal low-frequency accelerations of the solenoids. For measuring the tidal 
accelerations of the order of 10~ 13 cm/s 2 in a satellite the gravity gradiometer 
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[||] was proposed. However, the tidal acceleration of test bodies caused by 
low-frequency bursts is of the order of 1CT 16 -j- 1CT 18 cm/ s . 

2 A Peculiarity of Magnetically Interacting 
Superconducting Solenoids 

Consider a pair of the superconducting solenoids Pi and P2 in line, with 
points %\ and xi as the centers in the weightless state . If the solenoids carry 
the persistent currents I\ and I2 , their magnetic energies are U\ = l\Q\j2 
and U2 = I2Q2/2, where Qi and Q2 are magnetic fluxes in the solenoids. 
Since Qi = L\I\ + MI 2 and Q 2 = L 2 h + Ml\ , where L\ and L 2 are the 
inductances and M is the mutual inductance of the solenoids, the currents 
are 

h = (QiL 2 - Q 2 M)/D, I 2 = {Q 2 Lx - QxM)/D, (1) 

where D = LiL 2 — M 2 .The fluxes Q\ and Q2 are constants and M is a func- 
tion of the distance x = \x\ — x 2 \ between the solenoids centers. Therefore, 
l\ and I2 are the function of x and, when one of the conditions 

M = UQ2IQ1, M = L 2 Q 1 /Q 2 (2) 

is satisfied, the currents change their signs. It follows from eqs.(|l]) that the 
energy U = U\ + U2 of the system is given by 

U = (L 2 (Qi) 2 - 2MQ 1 Q 2 - L 1 (Q 2 ) 2 )/(2D). (3) 

The energy of the solenoids magnetic interaction is 

W = U - LxQl/2 - L 2 Q 2 2 /2. 

The Ampere force F = —dW/dx affects the solenoids. It follows from the 
eq.(|3|) that F = IiI 2 dM/dx , where I\ and I 2 are defined by eqs. ([I]) • 

The basic peculiarity of the interaction of superconducting circulating 
currents (The Kozorez effect ) is that the function W(x) has the minimum 
and the Ampere force is reduced to zero at the certain distance x = x , where 
one of the conditions ([TJ) is satisfied . 
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Indeed, F = U' M M' X ,where U' M = dU/dM and M' x = dMdx . However, 
the function U' M = [QiQ 2 M 2 + {L 2 Q\ + LiQ|)M - Q 1 Q 2 L 1 L 2 ]/ D 2 is equal 
to zero if one of the conditions (0) is satisfied. 

We have also Eft = U'i IM (M' x f + U' M M XX , where Eft = 0E£/as and 
ikft = dM x /dx .The mutual inductance M(x) is a monotonically decreasing 
function with an increasing distance x , and M'(x) 7^ 0. Let Q 2 < Qi. 
Since M < L\ , the first equation in (Q) is satisfied and under this condition 
U'mmM = {Qi) 2 /[(LiL 2 - M 2 )^} > .Then Eft(x ) > . If Q 2 > Q x , 
the second equation in (|2]) is satisfied and we again arrive at this conclusion. 
Therefore, the function W(x) has a minimum at the distance xq ■ If Q 2 < Qi 
, the energy W in the equilibrium position x = xq is W m i n = —Q^/i^L^. 
An illustrative example for the pair of the identical solenoids is given in Fig. 
1 . The parameters of the solenoids are : the inductances are l.lbHn. , the 
lengths are 5cm., the radiuses are 10cm., Qi = 1.15Wb ,Q 2 = Qi/25 . 
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Fig.l . The Force F(x) 
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3 The Magnetically Coupled Solenoids as a 
Detector of Tidal Accelerations 



Consider the properties of the system, formed by a pair of superconducting 
solenoids in weightless state, in the field of a gravitational wave with the 
frequency v extending during the time interval to perpendicularly to P1P2 
- direction. This system can be regarded as an oscillator, the small oscilla- 
tions of which in P1P2 - direction are described by the nonlinear differential 
equation 

mq + R(q) + pq = f(t). (4) 

In eq.(|4]) q = x — xo is a small deviation from the equilibrium position , 
q = dq/dt , q = dq/dt , R(q) is the air resistance, p = U"(x ) is the stiffness 
of the "magnetic spring", f(t) = ma g sm(ujt) at < t < t Q and f(t) = at 
t > t , m is the mass of the system , u = 2nv and a g = u 2 hx /2 is the 
amplitude of the tidal acceleration, caused by the gravitational wave. 

If 1-type superconductors are used in the solenoids, the air resistance 
R{q) is the key cause of the oscillations damping in the given system. (This 
assertion is argued in Section 4). Such a detector can be named an ideal 
detector. For an ideal gas the function R(q) = —bq\q\, where b = apS/2 , a is 
the aerodynamic factor of the resistance , p is the air density in the device, 
S is the crosssection of the solenoids orthogonally to the P1P2 - direction. 
If the pressure is 10~ 10 Torr in the hydrogen atmosphere, the temperature 
T = 4.2K and S = 10cm , the magnitude of b reaches 10~ u gm/cm . 

Since R(q) and a g are small quantities , we shall use the Bogolubov- 
Krilov method for the analysis of nonlinear equation (Q). 

Suppose, lj is close to the resonant frequency u = {p/m) 1 / 2 . We shall 
seek an approximate solution of eq.(f|) in the form 

q(t) = A(t) cos[ujt + (5) 

where A and d are slowly varied functions of time t. Let us replace R(q) in 
eq. (H) by the linear function X e q , where A minimises the function 

/(A) = f T [R(q) - X e q]dt (6) 
Jo 

(T = 27r). It is attained at A e = abAu, where a = 8/(37r). Denote X e /m by 
(3 . Then eq. (|J) takes the form of a linear differential equation 

q + (3q + {oJofq = a g sin(ujt), (7) 
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where, however, (3 is a function of A. 

It follows from eq.(3) that at Q\ > Q2 the stiffness p is given by 



p=[Q 1 M' x (x )} 2 /(L l D) (8) 

The typical values of uj are much less than 1Hz. Substituting eq. (||) into 
(£D, ignoring the terms A and (3 A and setting the coefficients of cos(c<jt) and 
sm(ut) equal to zero, we obtain the differential equations system for A and 
d : 

• a(3uj 2 A a , „ 

A + — + -^-cos(tf) = 9 

2m 2cu 

# ^ 2 -^o _ sinW = o (10) 
2u;m 2Am v ; v ; 

Setting A = and $ = and eliminating the phase $ , we obtain the 
equation 

A[{u 2 - uj 2 ) 2 + ab 2 uj 4 m- 2 A 2 } = [u A hx Q )/A (11) 

This equation gives the implicit function A = A(u) at stationary oscilla- 
tions. At the resonance (uj = u ) the amplitude is given by 

A = [mhx /(2ab)] 1/2 . (12) 

At uj = ujq and zero initial condition there exists solution A = A(t) with a 
constant phase $ = n : 

A(t) = A^tanh^Ao), (13) 

where To = c<j~ 1 [8m/(a6/ixo] 1 ^ 2 is the detector relaxation time .If in eq. (|l2l) 
m = 10 4 gm , x = 50cm , b = 10~ 3 gm/cm , h = 10~ 22 and uj = OAHz we find 
that the stationary amplitude of the detector response A m = 1.5 * 10~ 7 cm. 
However , it is impossible to observe such high amplitudes since relaxation 
time of the ideal detector ( tq is proportional to hr 1 ! 2 ) is too high. However, 
Ait) reaches 2 ■ 10~ 15 cm. already at the observation time t = 2.6 ■ 10 6 s (one 
month). 

At t <C to we obtain A(t) = (A m /To)t. Consequently, if the interval of a 
resonant gravitational-wave burst is T = 27r/a;, the detector response is 

q(t) = (nhx /2) cos(uj t) (14) 
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at t > T. Thus, the gravitational wave resonant burst gives rise to the long - 
duration , poorly damped oscillations with the frequency uo and the amplitude 
7r/iXo/2 . At h — 10~ 17 and Xq = 50cm the amplitude is 7.8 • 10~ ie cm. 

If the resonant frequency u Q is much less than u , then a numerical solu- 
tion of eq.fl^) shows that the response q(t) reaches 1.7* 10~ 15 cm at t — 10s 
and further on slightly varies in time . 



4 Measurement noises 

The detector under study is a nonlinear oscillator. Let us find the variance 
of its thermal fluctuations that is other than the one of a linear oscillator. 
Consider first the linear oscillator describred by the differential equation of 
the form q + Pq + u^q = f(t). Suppose, the oscillator relaxation time is much 
longer than the measurements time, i.e. fit <C 1. In that case , acording to a 
rigorous solution by Chandrasekhar || the position q variance (u q ) 2 is given 
by 



sin(2u;ot) 
2uj t 

and the velocity v = q variance is given by 



2 kT n 

a\ = ^(3t 



+ o((3t 2 ) (15) 



2 kT n 

a 2 v = —fit 
m 



1 + Sin(2 " ot) + o{^) 
2ujt 



(16) 



The expressions between the brackets are essential at low resonant fre- 
quencies of the oscillator. 

In the case under consideration for an approximate analysis of the ther- 
mal fluctuations it is expedient to take the statistic linearization of the func- 
tion R(q). Namely, we put R(q) = Xi e q , where Ai e minimises the function 
A(Aie) = ([R(q) — Ai e <?]) • ( The symbol ( ) denotes an ensemble mean). 
Now the random oscillations of the detector, caused by the thermal noise, 
are given approximately by the differential equation q + bq + u 2 q = £(£), 
where (3 = \\ e jvn and is a Gaussian noise. 

The desired value of A ie is Ai e =< qR(q) > / < q 2 > . 

Assuming that the distribution function is approximately Gaussian, we 
find that (3 = 4(27r) -1 / 2 m~ 1 a v . Thus, (3 is the function of a v . 

Now we are considering eq. ( |T6"D as an equality whence the function a v (t) 
can be found. Substituting the above expression of (3 into eq.([T6|) and found 
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function a v (t) - into eq.(R) we obtain, finally, that the value of the variance 

2 . 



Ab 2 kT 



sin (2uj t) 



2u 2 t) 2 



(17) 



For example, if m = 10 4 gm , uq = 0.01s, t = 100s , b = 10~ A gm/cm, 
then the root-mean-square magnitude of the thermal fluctuations is : < 
q 2 > l l 2 = a q = 10~ 20 cm. Thus, the detector has a very low level of the 
thermal noise, that is much less than the expected detector response to the 
gravitational wave bursts. 

The inhomogeneity of the Earth and spacecraft gravitational fields lead 
to more serious problems. The difference in the solenoids gravitati- onal 
accelerations approximately equal to x , where U^. = d 2 U / dx l dx k and 
U is the Earth gravitational potential. The variations of JJ^. during the 
orbital motion of the spacecraft are well beyond the tidal accelerations a g 
in eq.(4). This problem can be solved by choosing a geostationary or very 
distant from the Earth spacecraft orbit. 

It is necessary to take into account the fluctuations in the gravity gradient 
within the spacecraft too. For example, the 1cm shift in the position of the 
5gm mass at its distance 2m from the solenoids along the -P1-P2 -direction 
causes the variation 4-10 _14 cm/s 2 in the relative acceleration of the solenoids. 

A little part of the inductance (the motional inductance) is a function 
of the temperature T. So, variations in T cause the variations 6L in the 
inductance L. (In the above illustrative example at the wire radius 0.05cm 
dL/L is about 10~ 16 at T = 2K and 10~ 20 at T = 0.1K ). 

For definiteness, assume that Q% < Q2 ■ Then under the equilibrium 
position Q2L1 —QiM = and D w L\L 2 ■ If at moment t = a temperature 
variation begins, then, according to eqs. (2), at t > we have I\ w 7° + 
I^(5L 2 /L 2 ) and I 2 w I^SLt/L^ ,where If = Q 1 /L l and 1° = Q 2 /L 2 . Now 
the force of the interaction between the solenoids in absence of exterior forces 
is other than zero : 

F T = I° 1 I^SL 1 /L 1 )M , (x ), (18) 

whereM'(xo) = dM/dx . The force Ft gives rise to slow oscillations of the 
solenoids relative to the equilibrium position xo ■ 

Denote the acceleration Ft/tti by ar ■ Obviously, it is necessary to select 
the detector parameters so that during the measurements time the inequality 
a g < ax is valid. We do not analyze the problem of optimization of the 
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detector parameters here. However, it should be noted that it is easiest of all 
to obey the above inequality if the detector resonant frequency u is much 
less than the frequency uo of the wave . This happens to be the case for 
the considered above illustrative example (Fig.l) at Qi = 1.15 • 10 — 1 In 
that case u = 2.6 • lO^s" 1 , 7° = 0.1 A and J 2 ° = 7°/25 ). At T = 2K 
and 5L/L = 10~ 16 the acceleration = 5 • l0~ ls cm/ s . Meanwhile, at 
v = 0.1Hz and h = 10~ 17 the gravitational acceleration a g = 4 • 10~ 17 cm/s. 
(The detector response to short bursts found by the the numerical solution 
of eq.(4) is about 7 • lCT 16 cm at to = 10s). 

It follows from eq.(8) that the stiffness p = [I^M' (xq)] 2 / L 2 . Hence, the 
temperature fluctuations cause the resonant frequency variations 5ujq . They 
are given by 5ujq/uiq = —SLi/ (2Li). These variations may be essential at the 
resonant detection of a continuous signal. 

Under oscillations of the solenoids relative to the equilibrium position 
the persistent currents in the solenoids are not constants ( I = dl/dt ^ ). 
From the viewpoint of the two-fluid model the intensity of the supercurrent 
is J = en s v s , where e is the charge of the electron, and n s and v s are the 
concentration and velocity of the superconducting electrons , correspondingly. 
As a result of the acceleration of the superconducting electrons an electric 
field E appears inside the superconducor that can be found from the equality 
mv s = —eE , where v s = dv s /dt. If the time -depending current is of the 
form Jexp(— iujt) , then E = —[imuJ/(e 2 n s )}. This electric field causes a 
motion of the normal electrons according to the following equation 

md(v n )/dt+ {m/r )(v n ) = -eE, 

where (v n ) is the mean velocity of the normal electrons and r is the relax- 
ations time (usually about 10 _13 s). Then, there is also a normal current in 
the solenoids with the intensity J n = —en n (v n ) , where n n is the normal 
electrons concentration. Solving the above equation of motion we find that 
at low frequencies uj the normal current is given by 

In = {n n /n s )uJTl° a sva(ujt). (19) 

Let Qi < Q2 ■ The amplitude I™ ax of the supercurrent in the solenoid 
P2 caused by gravitational-wave bursts is equal to I' 2 A , where A is the 
amplitude of the detector response and I' 2 = dl 2 /dq . It follows from eq.(l) 
that I[ = (1%/L^M' and V 2 = (I°/L 2 )M' .The typical value of V 2 is of the 
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order of 10~ 2 -j- 10~ 4 Aj cm. Since A is about 10 _16 cm, the typical value of 
jmax p roves ^ Q about 10~ 19 A . Setting in (19) n n /n s = 1/5 and v = 
0.1Hz we find that the amplitude of the normal current is 10~ 34 A During 
time dt the energy dissipation is RI 2 dt, where R is the superconductor wire 
normal resistance . If R = 1Q , then during the time T = the energy 

dissipation is about 1CT 67 J . At the same time the tidal forces work is about 
ma g A = 10~ 36 J . 

Consider another question: ought we to include in eq.(10) in addition to 
the term R(q) an effective force F e f , which describes the energy dissipation 
caused by the normal current ? In principle, it can be done since according to 
eq.(19) I n is of the form (n n /n s )rl s (t). Hence, RI 2 dt can be written as F e f<jdt 
, where F e f = Rr 2 (dl /dq) 2 q . Let us compare F e f with the force R(q). The 
value of q is approximately ujA . At low frequencies and at b ^> 10~ 10 gm/cm 
the following inequality is fulfilled: R(q) < F ef , since Rt 2 (I') 2 <ti bAu . 
(Because of too longe relaxation time r , small values of b ought not to be 
used ). 

Thus, at very small amplitude oscillations and low frequencies the energy 
dissipation caused by the normal current is insignificant. 

Consider briefly the physical principles of the solenoids relative shift mea- 
surement caused by the gravitational-wave bursts. For definiteness, assume 
that Q2 < Qi ■ At the equilibrium position the supercurrent I2 = 0. The de- 
viations in the equilibrium position cause the change in the proper magnetic 
flux of the solenoid P 2 is SQ = 5I 2 L 2) where 5I 2 = I' 2 q ■ Thus, the value 5Q 2 
is given by 

5Q 2 = l\M'q (20) 

Suppose that the distance between the near solenoids ends is less than 
lcm at the distance between the solenoids centers of 50cm . In that case the 
value M' reaches 10~ 2 Gn/cm or more than that. At I® = 0.1A dQ 2 is about 
10~ 3 $ > where $ = 2lO~ 15 Wb is the magnetic flux quantum. 

A method of dQ 2 measuring is shown in Fig. 2 . In this schematic dia- 
gram S is a superconducting quantum- interferometer device (SQUID) that 
is attached to the solenoid P 2 and coupled with P 2 inductively by a flux 
transformer T. 
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Pi P2 




Fig.2 

The minimal magnetic flux measured by the SQUID with the T is given 
by |§ 5Qmin = 2[2L a £(z/)] ' /iV a , where L a is the antenna inductance , N a 
is its turns number and E{y) is the SQUID noise expressed as an input en- 
ergy resolving at the frequency v . At the frequencies v < 0.1 Hz the value 
£{y) = (lO-^/v^J/iHz)- 1 / 2 . At v = 0.1Hz 5Q min = 10- 4 $ or less than 
that. Thus, 5Q min < 5Q 2 ■ 

We mean that the solenoids Pi and P 2 are inside a supeconducting shield. 
The results of the SQUID measurements are transmitted by radio by means of 
a conversion "voltage - frequency" and by using an isotropic active antenna. 
Such a method of the solenoids shift measuring is insensitive to micromete- 
orites impacts and other forces affecting the spacecraft. 

In the SQUID circuit there exists noise magnetic flux $ ra of the order 
of 10 _5 $o or less than that. Because of the inductive coupling between the 
SQUID and the solenoid P2 a noise current appears in the latter that is less 
than $/L 2 = 10 -20 v4 . It is much less than the noise effect upon the detector 
is a small value. 

A number of the other disturbing forces in the space laboratory is con- 
sidered in [11 ] . The author would like to thank C.W.F. Everitt and 
P.W.Worden Jr. for the materials on the STEP. 



10 



References 

[1] K.S. Thome , Rev. Mod. Phys. 52 , 285 (1980) 

[2] V.B.Braginsky ., K.S.Thorne . Nature 316 , 610 (1985) 

[3] H.A.Chan and H.J.Paik ,. Phys. Rev. D 35, 3551 (1987). 

[4] V.V.Kozorez , Dynamics of the magnetic interacting bodies systems 
Kiev (1981) (In Russian). 

[5] A.Blaquier , Nonlinear systems analysis, New York and London (1966). 

[6] S.Chadrasekhar , Stochastic problems in Physics and Astronomy, 
Rev.Mod. Phys. 15 , 1 (1943). 

[7] T.Van Duser and C.W.Terner ,. Principles of superconductive devices 
and circuits , Elsvier North Holland, Inc. (1981). 

[8] Superconductor Application : SQUIDs and Machines. Ed. by 
B.B.Schwartz and S.Foner , Plenum Press, New York, (1977). 

[9] J.P.Blaser at al, STEP (Satellite Test of the Equivalence Principle), 
Report on the phase a study. ESA and NASA , SCI(93)4 March 1993. 



11 



